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We plan to determine the bulk viscosity of xenon 
10 times closer [in reduced temperature r = ( T-T c )/T c ] to 
its liquid-vapor critical point than ever before. ( T c is the 
critical temperature.) To do so, we must measure the 
dispersion and attenuation of sound at frequencies 1/100 
of those used previously. In general, sound attenuation 
has contributions from the bulk viscosity acting 
throughout the volume of the xenon as well as 
contributions from the thermal conductivity and the shear 
viscosity acting within thin thermoacoustic boundary 
layers at the interface between the xenon and the solid 
walls of the resonator. Thus, we can determine the bulk 
viscosity only when the boundary layer attenuation is 
small and well understood. We present a comparison of 
calculations and measurements of sound attenuation in 
the acoustic boundary layer of xenon near its liquid-vapor 
critical point. 

We used a novel, compact, acoustic resonator 
designed specifically for these measurements (Fig. 1). We 
measured the frequency response of this resonator filled 
with xenon at its critical density p c in the reduced 
temperature range 10 3 < T < 10" 1 . From the frequency- 
response data, we obtained the resonance frequency and 
the attenuation for six resonant modes in the range 
0.10 kHz </< 7.5 kHz. Using the known thermo-physical 
properties of xenon, we predict that the attenuation at the 
boundary first increases and then saturates when the 
effusivity of the xenon exceeds that of the solid. [The ' ^ 

effusivity is s = (pCpX T ) , where Cp is the isobaric Figure 1. Acoustic resonator 

specific heat and At is the thermal conductivity.] The 

model correctly predicts (±1.0%) the quality factors Q of resonances measured in a steel 
resonator (£ ss = 6400 kg-K' s’ ); it also predicts the observed increase of the Q, by up to a 
factor of 8, when the resonator is coated with a polymer (s pr = 370 kg-K' 1 -s' 5 2 ). The measured 
acoustic dissipation in near-critical xenon shows a prominent plateau for the Helmholtz mode 
and, to a lesser extent, for the longitudinal modes (Fig. 2). These results are the first direct 
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evidence of thermal boundary dissipation being limited by the thermophysical properties of the 
solid wall. 





Figure 2. (a) Dissipation due to bulk viscosity 

and the thermal boundary layer versus reduced 
temperature. Theory for thermal boundary 

dissipation ( -). The fractional deviations of 

the measured dissipation from theory versus 
reduced temperature for all the modes in the steel 
(b) and polymer-coated (c) resonators. Only data 
for which the dissipation from bulk viscosity was 
less than 0.8 % of the total dissipation is shown. 


We observed an 8-fold reduction in the 
thermal boundary dissipation (as predicted) 
after coating a steel resonator with a low- 
effusivity polymer. (See Fig. 2a.) Such a 
reduction in thermal dissipation will be 
necessary in order to measure the low- 
frequency dissipation from bulk viscosity close 
to the critical point in future experiments. 

To analyze the data, we formulated a theory 
for acoustic dissipation in a fluid that is 
bounded by a rigid wall with the condition that 
the temperature and heat flux across the 
boundary be continuous. The theory is valid 
both near to and far from the critical point; it 
includes volume and surface dissipation from 
thermal conduction, shear viscosity, and bulk 
viscosity. 

We measured the speed of sound in xenon 
as a function of reduced temperature in the 
range 0.0006 < t < 0.03. Using the measured 
sound speed and attenuation, the isothermal 
susceptibility from Ref. [1], and ( dP/dT) p from 
Ref. [2], we derived Cy, Cp, and the 
background terms for the thermal conductivity. 
These derived properties were used to 
determine the amplitude A + = 18.01 for the 
singular part of Cy (with a = 0.110) and the 
correlation length amplitude £ 0 + = 0.1866 nm. 

Although the data do not extend far into the 
asymptotic region, the agreement with other 
values is remarkable. 

The onset of bulk viscosity is evident from 
the excess dissipation (over the thermal and 
viscous dissipation) indicated by the upturn at 


low r of the measured dissipation shown in 
Fig. 2a. It is also evident that the dissipation 
from the bulk viscosity increases as the frequency increases in qualitative agreement with theory. 
The onset of bulk viscosity is more evident in the polymer-coated resonator because the thermal 
dissipation is smaller than it is in the bare-steel resonator. 

[1] H. Giittinger and D.S. Canned, Phys. Rev. A24, 3188-3201 (1981). 

[2] H.L. Swinney and D.L. Henry, Phys Rev A 8 , 2566-2617 (1973). 


NAS A/CP— 2004-2 1 3205/VOL2 


101 





NAS A/CP— 2004-2 1 3205/VOL2 


Progress on acoustic measurements of 
the bulk viscosity of near-critical xenon (BVX) 

Keith A. Gillis (Co-1), Iosif I. Shinder, Michael R. Moldover (PI) 

National Institute of Standards and Technology, Gaithersburg, MD 20899 

and Gregory A. Zimmerli (Co-1) 

NASA Glenn Research Center, Cleveland, OH 44135 

We plan to determine the bulk viscosity of xenon 10 times closer [in reduced temperature 
t- (T-T c )/T c ] to its liquid-vapor critical point than ever before. ( T c is the critical temperature.) 
To do so, we must measure the dispersion and attenuation of sound at frequencies 1/100 of 
those used previously. We will determine the bulk viscosity C, (also “second viscosity” or 
“dilational viscosity”) from measurements of sound attenuation a A and dispersion using a 
novel, compact, acoustic resonator designed specifically for this purpose. 

In general, sound attenuation has contributions from the bulk viscosity acting throughout the 
volume of the xenon and contributions from the thermal conductivity and the shear viscosity 
acting within thin boundary layers near the resonator wall. Thus, we can best determine the 
bulk viscosity when the boundary layer attenuation is small and well understood. We present 
a comparison of calculations and measurements of sound attenuation in the acoustic 
boundary layer of xenon near its liquid-vapor critical point. 
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The liquid-vapor critical point 

The critical point is the highest temperature and pressure at which liquid and 
vapor can coexist in the same container. For xenon T c =16.584 °C, P c = 5.84 MPa 



Temperature (°C) 

As the critical point is approached: the sound speed — > 0; the bulk viscosity £ diverges; y= C t JC v diverges. 
Attenuation in the thermal boundary layer is proportional to y and therefore poses a potential problem. 
Conventional acoustic theory assumes y~ 1 and attenuation is a small perturbation. 
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Accomplishments 

• We extended acoustic theory to include near-critical fluids, 
heat flow in solid wall, and bulk viscosity. 

• We were the first to observe direct evidence of thermal 
boundary dissipation being limited by thermophysical 
properties of the solid wall. 

• We demonstrated an 8-fold reduction in thermal boundary 
dissipation after coating steel with a polymer. 

• We derived C v , C p , and background thermal conductivity 
from our measured speed of sound. 

• We observed the onset of dissipation due to bulk viscosity. 
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Helmholtz mode 


O 

Ltl 



cot 



k H = 0.1226cm" 1 => X = 51.0cm ~ 34 


In it’s fundamental acoustic mode, the gas oscillates between two chambers through a small tube (radius = r d and 
length = LJ. The wavelength of this mode is about 10 times the size of the resonator or 34 times the tube length. 
We measure the response as a function of frequency to determine the resonance frequency and Q. For xenon at 
its critical density, the frequency of this mode varies from 127 Hz at 3 mK above T c to 240 Hz at 3 K above T c . 


We use resonators of this type (called double-Helmholtz resonators) to measure the viscosity of dilute gases. 
See K.A.Gillis, J.B. Mehl, and M.R. Moldover, “Theory of the Greenspan Viscometer,” JASA 114, 166-173 (2003) 
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Longitudinal modes 



COt cot + 71 

k Ll = 0.6536cm' 1 => A = 9.6cm = 2 L c 

We deduce the sound speed and the acoustic attenuation from the resonance frequency and Q for a given mode. 

To measure dispersion, we also study the first 5 longitudinal modes of the longer chamber. Like the Helmholtz 
mode these modes are non-degenerate and well separated because of the asymmetric resonator design. All 6 
modes span a factor of 27 in frequency at a given temperature. The acoustic data span a factor of 60 in frequency 
over the full temperature range. 
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The 2 chambers and the connecting tube were machined as one piece from a block of stainless steel. The high 
rigidity is crucial for studies of dense fluids. Piezoceramic disks were attached to the outside of the flanges in 
order to drive and detect sound waves in the fluid. Two similar resonators were used. One resonator was coated 
with a polymer to reduce the thermal conduction between the xenon and the resonator wall. The other resonator 
was left as bare steel. 
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Acoustical Spectrum of Resonator 
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The graph shows the amplitude of the detector signal as a function of frequency. The Helmholtz mode and the 
longitudinal modes of the chamber (L1-L5) are clearly visible and well-separated. Isolated, non-degenerate 
modes are necessary for accurate modeling. 
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From dilute gas to the critical isochore 


O 

VO 


0.1 


g 


o.oi 


o.ooi 



i 


10 100 1000 10000 


density (kg/m 3 ) 


This graph shows the measured dissipation (1/Q) as the density increases from ambient conditions to the critical 
density (1116 kg/m 3 ) at constant temperature (blue curve). The contributions from viscosity and thermal 
conductivity within the boundary layers, calculated from the model, are shown in green and red. Viscous 
dissipation dominates at low density. Thermal dissipation dominates at high density due to the large y. The black 
dashed line shows the dramatic increase in 1/Q as the temperature is lowered to T c along the critical isochore. 
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Dependence of dissipation on thermal properties of the wall 
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To lowest order, the dissipation is due to the sum of 3 contributions: thermal and viscous dissipation at 
the interface between xenon and the wall plus a volume dissipation due to bulk viscosity. The thermal 
dissipation reaches its maximum value when the ratio 1 } of the effusivity of xenon e to the effusivity of the 
solid wall e sec equals one. Lower £ soljd means the thermal dissipation saturates at a lower value farther 
from T c . See K.A. Gillis, l.l. Shinder, and M.R. Moldover, “Thermoacoustic boundary layers near the liquid- 
vapor critical point,” (in press Phys. Rev. E, July 2004). 
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Theoretical predictions 
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Predicted dissipation as a function of reduced temperature in bare steel and polymer-coated resonators. The 
crosses indicate where the effusivity e of xenon equals the effusivity of the solid. The surface dissipation from 
viscosity and heat conduction is shown with short dashed curves. The contribution from bulk viscosity £ in the 
volume is shown with a long dashed curve. The relaxation of fluctuations gives rise to bulk viscosity and has 
characteristic time r^. The condition yco£ I pc 2 ~1 indicates where the boundary layer becomes highly compressible. 
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Comparison of measurements with theory 

Xenon near its critical point 
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(a) Dissipation due to bulk viscosity and the thermal boundary layer versus reduced temperature. 

Theory for thermal boundary dissipation ( ). The fractional deviations of the measured 

dissipation from theory versus reduced temperature for all the modes in the steel (b) and polymer- 
coated (c) resonators. Only data for which the dissipation from bulk viscosity was less than 0.8% of 
the total dissipation is shown. 
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Speed of sound in xenon 
near its critical point 


(a) The speed of sound c HS determined 
from the Helmholtz mode with bare 
stainless steel plotted versus the 
reduced temperature. The speed of 
sound measured by Kline and Carome 
at 6 kHz. 


(b) Deviations of the speed of sound 
determined from each mode from c HS 
shown in the upper graph. 


Consistency among modes -0.01% 
Consistency between resonators -0.1% 


